Diquark Bose Condensates in High Density Matter and Instantons 
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Instantons lead to strong correlations between up and down quarks with spin zero and anti- 
symmetric color wave functions. In cold and dense matter, Ub > Uc — lfm~'^ and T < Tc ~ 50 
MeV, these pairs Bose-condense, replacing the usual {qq) condensate and restoring chiral symmetry. 
At high density, the ground state is a color superconductor in which diquarks play the role of Cooper 
pairs. An interesting toy model is provided by QCD with two colors: it has a particle-anti-particle 
symmetry which relates (qq) and (qq) condensates. 
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The properties of hadronic matter under extreme con- 
ditions are subject to intense theoretical studies, nu- 
merical simulations on the lattice, and experimental ef- 
forts using high energy heavy-ion collisions. Substantial 
progress has been made with respect to high tempera- 
ture QCD matter, while the (more difficult) problem of 
cold dense matter is much less understood. This is partly 
due to the fact that up to now, lattice simulations have 
not been able to overcome the numerical problems associ- 
ated with the fact that at non-zero chemical potential the 
fermionic determinant is complex (see ^ for a recent re- 
view). This is unfortunate, because the phase structure 
of dense matter may be very rich. Several intermedi- 
ate phases between nuclear and quark matter have been 
proposed, e.g. states with mesonic (pion or kaon |^) 
condensates. 

Due to the asymptotic freedom of QCD one expects 
very dense matter to resemble an ideal Fermi gas of 
quarks. The system is quite similar to a cool electron 
plasma, with Debye-screening of color fields at momen- 
tum scales p < Mo ~ collective plasmon excitations, 
etc. Since the Coulomb interaction between quarks of 
different colors is attractive, it was realized early on that 
the quark plasma should be a superconductor, due to the 
formation of Cooper pairs on the Fermi surface Q. The 
magnitude of the corresponding gap A and the critical 
temperature Tc were estimated to be in the MeV range. 

In this article we show that non-perturbative effects 
lead to diquark condensates with A,Tc about two or- 
ders of magnitude larger. These condensates are gener- 
ated by the instanton-induced interactions between light 
quarks For two flavors (up and down) the (qq) inter- 
action is 
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where we have added the interaction in the direct and 
exchange channels and dropped color octet terms. Nc 
is the number of colors and — (r, i) is an isospin 
matrix. We will specify the coupling constant G below. 
There is pervasive evidence for the importance of this 
interaction from (i) phenomenological studies of current 
correlation functions in QCD, (ii) the success of hadronic 
spectroscopy in the instanton liquid model, and (lit) stud- 



ies of instantons and their effects on the lattice, see ^ 
for a review of these issues. 

The result (|l]) can be Fierz-rearranged into a (qq) in- 
teraction. One obtains 



C = G 



+ (^^CT2Al75V')(^T2Al75CV5^)] 

+ 32jV,(A^, + 1) 

Here, C is the charge conjugation matrix, T2 is the anti- 
symmetric Pauli matrix, Xa,s a-re the anti-symmetric 
(color 3) and symmetric (color 6) color generators. The 
effective lagrangian (^) provides a strong attractive in- 
teraction between an up and a down quark with anti- 
parallel spins {J^ = 0+) in the color anti-triplet channel, 
and a repulsive interaction in the 0~ channel. 0"*" quark 
pairs couple to the diquark current S^^ = eabcU^C^^dc- 
Phenomenological implications of the instanton-induced 
interaction in this channel were first discussed in con- 
nection with spin-dependent forces in baryons [Q, chal- 
lenging the conventional wisdom that spin splittings are 
due to one-gluon exchange. The importance of diquark 
degrees of freedom is perhaps most obvious in baryons 
that contain one very heavy quark. Quantitative studies 
of instanton effects in baryon spectroscopy (both light 
and heavy- light systems) were done in ||]. The conclu- 
sion was that the instanton effects encoded in (H) are 
indeed strong enough to reproduce the observed spin- 
splittings, and that the nucleon has a very large over- 
lap with the current eabc{u^C"f^d'^)u'^ = S^gU""- Since 
there is no confinement in the instanton model, one can 
compare the diquark mass to the two constituent quark 
threshold. The result is a deeply bound scalar diquark 
2mq — msdq — 200 — 300 MeV , whereas all other chan- 
nels (vectors and axial- vectors, color 6 diquarks, etc.) are 
at most weakly bound. 

The possible role of diquark clusters in quark matter 
was discussed in . It was noted that a loosely bound 
"third" quark in the nucleon may find a partner in dense 
matter. However, as we show below, this effect is less im- 
portant than Bose condensation. In general, even if we 
focus only on the instanton-induced interaction, the situ- 
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ation is quite involved due to the competing attraction in 
the qq channel. At small density, the attraction in the qq 
channel is stronger than the one in the qq channel (com- 
pare (^) with (^)), qq pairs condense and chiral symmetry 
is broken. At large density, Pauli-blocking suppresses the 
qq interaction, but the qq interaction leads to an instabil- 
ity near the Fermi surface. The mechanisms for qq and 
qq condensation are quite similar, both are based on the 
formation of a gap in the fermion spectrum at the surface 
of the Dirac/Fermi sea, respectively. We will defer an at- 
tack on the coupled problem to a separate publication, 
and in this work only consider the low and high density 
parts separately. 

At low baryon density, we ignore modifications of the 
vacuum and keep (qq) fixed. In this limit instanton-based 
models realize the minimum of the energy per quark in 
scalar diquarks, not in nucleons. Only by taking confine- 
ment into account, one finds that at low density the or- 
dinary nuclear matter is favored over a diquark Bose gas. 
At high densities we ignore interactions with anti-quarks, 
assume a deconfined chirally symmetric quark matter, 
and calculate how instantons operate at the Fermi sur- 
face. For pedagogical reasons, let us first discuss the 
special case of QCD with two colors Nc — 2. 

2. In the massless case Nc=2 QCD has an addi- 
tional Pauli-Giirsey symmetry (PGSY) [ p]Jl^ , which 
mixes quarks with anti-quarks. As a result, diquarks 
(=baryons) are color singlet states that are degener- 
ate with the corresponding mesons. Chiral symmetry 
breaking then implies that some diquarks are Goldstone 
bosons, that means their mass vanishes in the chiral limit 
{niq = 0). Their properties follow directly from the 
PGSY. The total number of Goldstone modes is 13-1^] 
2iV| — Nf — 1. For Nf=l there is no Goldstone boson 
(the U{1) symmetry is anomalous). For Nc = Nf = 2 
one finds five: 3 pions, the scalar diquark S and its 
anti-particle S. There is a nice continuity in going from 
iVc = 2 to 3: the scalar diquark goes from being massless 
to a deeply bound state! 

The coset of the full group over the unbroken one is 
K = 5L/(4)/5p(4) = SOi6)/SO{5) = S^, i.e. the effec- 
tive chiral Lagrangian describing the Goldstone modes is 
a sigma model, in which the usual "chiral circle" is re- 
placed by a 5-dimcnsional sphere. The ordinary vacuum 
of the massless Nc — Nf = 2 theory with a non-zero 
"mesonic" condensate (qq) can simply be rotated into 
states with finite diquark condensate. This rotation costs 



no energy, so the chemical potential fi = dE/duq 
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At the point where the {qq) condensate is as large as the 
original {qq) condensate, chiral symmetry is restored, but 
there are still 5 Goldstone modes, which are now 3 pions, 
sigma and S. 

For non-zero quark mass the qualitative picture can be 
understood using the corresponding linear sigma model. 
The potential 

V = \ (TT^+a^ + S^ + S^-w^)^ -Aa-^? {S'^ + S'^) (3) 



includes the diquark chemical potential /i and the chirally 
asymmetric mass term A. At /i = the Goldstone masses 
are — Ajv^ and — 8Aw^. For non-zero /i we 
can determine the {qq) and {qq) condensates (ct) and {S) 
using the mean field approximation. We find 
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Below the critical chemical potential [ic — Tfigl ^f2^ {a) is 
constant and {S) = 0. Above /ic, (5*) increases as 
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and (tr) = m'^v/{2^^). The energy density is e 
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3^2/(4^)2, compared to e = — m^t;^ -I- mg/(16A) for the 
normal vacuum. 

Unlike real QCD, Nc = 2 gauge theory is straightfor- 
ward to simulate on the lattice, since the fermion deter- 
minant remains real for /i 7^ 0. With the exception of 
some early work using small lattices and the strong cou- 
pling expansion ||l^, few studies have taken advantage 
of this. Numerical studies of the instanton model for 
Nc = 2 at finite density |17| are consistent with the sce- 



nario found above, at large density the {qq) condensate 
is replaced by a {qq) condensate. 

3. Let us return to the low density limit of real QCD 
with three colors, ignoring strangeness. For definiteness, 
we consider neutron (udd) matter relevant for stars [p^ in 
which ud diquarks and d quarks compensate each others 
color and electric charges. Since scalar diquarks are color 
anti-triplets, a Bose condensate will select a direction in 
color space. If we label this direction red (fc=3), our ud 
diquarks are made of blue and green (fc=l,2) quarks only. 
Their properties are strongly modified by the diquark 
condensate, while the third quark-type (red d) is basically 
unaffected. 

Diquarks are favored over a quark Fermi gas due to 
both their binding mg < 2mcff as well as their Bose 
character. However, it is erroneous to conclude that an 
infinite number of diquarks condenses in the p = state: 
diquarks are composite objects and, like nucleons, their 
interaction should have a repulsive core. We account 
for this by introducing a scattering length |19| a ~ 0.3 
fm into the expression for the energy per quark in the 
diquark Bose gas: 



n„ m 
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where ns is the S diquark density. The first term is the 
mean field result, and the second term comes from non- 
condensed diquarks |po| . The repulsion makes the pure 
diquark gas less favorable than the optimal Bose-Fermi 
mixture. Results of our calculations are shown in Fig.^, 
in which matter consists of (i) a Bose gas of S diquarks in 
chemical equilibrium with (ii) a Fermi gas of blue/green 
quarks, and neutralized in color and electric charge by an 
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appropriate amount of (iii) red quarks. For definiteness, 
we use the di-/quark masses of 500 MeV and 400 MeV, 
respectively. 

To account for confinement we add color strings to our 
picture. For light hadrons, lattice calculations suggest as 
potential of the form V{R) = K\R - i?o|6'(i? - Rq) with 
the usual string tension ~ 1 GeV/fm but smoothed 
out at the origin, here represented by Rq ~ 0.7 fm (see 
example in Q). We evaluate the average energy of a 

— 1/3 

string using this schematic potential with R = n^^^ 
where ritot is the total density of all quarks and diquarks. 
As can be seen from Fig.|l|, strings preserve nuclear mat- 
ter at low density. A diquark/quark/string mixture has a 
shallow (meta-stable) minimum leading to a mixed phase 
at n « 0.5-1 fm~^. It may be an artefact of the crude 
model: but most recent lattice results have some in- 
dications for such behaviour. Note a significant gain in 
total energy relative to quark matter (upper grey curve) . 
The critical temperature for Bose condensation can be 
roughly estimated from Einstein's ideal gas expression, 
Tc = 3.31n^^'^/m, which is about 100 MeV at the cross- 
ing point. 
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quarks explicitly, but parametrize the /x-dependence of 
(qq) and the constituent quark masses rriq by the factor 
(/i^ — M'^)/Mc foi' M < — 500 MeV (resembling a lin- 
ear decrease in density), where the precise value of the 
critical chemical potential for chiral restoration, fic, will 
be determined in future work. The gap equation in the 
most attractive channel mdjeij^ then reads 



1 



(2^) 



p dp 



tanh(ep(A)/2r) 



(7) 



with ep(A) = [(Wp-/l)2-HA(/i,r)2]2, ^2 ^p2_^i 



II — m^. For fj. > fic the Debye screening of the 
instanton fields becomes effective , which we account 
for in the coupling constant of the effective fermionic in- 
teraction: 

ge//(/i) = C (8^V.9')'/ dp p exp[-8^V5(p)'] 



X eicp[-Nfp^{p^ - Mc)^'(m - A^c)] exp[-Ap2 



(8) 



The last exponential factor is included to provide a cutoff 
at large p. Using the gap equation in the quark-antiquark 
channel we fit the value of the constant A in vacuum 
{T=p—Q) to a constituent quark mass mq ~ 400 MeV. 
In QCD with 3 flavors geff{p) for ud quarks is addition- 
ally reduced by the factor (m° -I- m*{p,))/{mP^ + m*(0)) 
due to the decreasing (ss)-related contribution to 
(m:(0)~200 MeV). Fig. 2 shows the gap A(/i,T = 0) 
with Nf=3 for two different values of the critical chemi- 
cal potential. At large p the gap is strongly suppressed 
by screening effects, while at small p it is reduced due to 
the decrease of the density of states at the Fermi surface 
(here the approach discussed in sect. 3 is more appropri- 
ate). The maximum gap is approximately linear in the 
critical density and reaches 50-100 MeV. 




FIG. 1. Energy per baryon (in MeV) versus baryon charge 
density rib (fm~"^). 

4. We now turn to the high density limit, with most 
of the quarks forming a Fermi gas and instanton-induced 
forces operating only near the Fermi surface. The width 
of this zone (the analogue of the Debye frequency in a 
phonon superconductor) is governed by the instanton 
form factor. The fermion zero modes that lead to the 
interaction (2) and determine the form factor depend on 
density They oscillate in the spatial direction, re- 

flecting the presence of a Fermi surface, and are long 
range in time direction. As a result, the diagrams in the 
qq channel are infrared divergent, making the Fermi sur- 
face unstable and leading to the formation of a gap. 
The full problem of evaluating the size of the gap is quite 
involved. Instead, we will consider a BCS-type gap equa- 
tion where the form factor is represented by a fixed cutoff 
A~G.3 GeV. We do not include interactions with anti- 
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n, [fm"'] 

FIG. 2. The gap A(/i,r=0) for pc=OA, 0.5 GeV (dashed 
and full line, respectively) and critical temperature Tc (dash- 
dotted line, ^c=0.5 GeV) versus baryon charge density ni,. 

Also plotted is the critical temperature Tc at which su- 
perconductivity disappears. The largest Tc occurs at 
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slightly higher densities than the largest gap, but at fixed 
chemical potential the BCS relation Aq = 7rTc/7 approx- 
imately holds (where 7 ~ 1.78 is Euler's constant). 

With fi exceeding the strange quark mass of 
rrig ~150 MeV, s-quarks are also present, which opens the 
possibility of forming new condensates. The instanton- 
induced pairing coupling constant in the ud channel is 
©(mO), while it is ©(mO) {0{ml)) for us {ds). There- 
fore strange-light quark pairing is strongly suppressed, 
presumably superseded by perturbative Coulomb inter- 
actions 0], which leads to A~l MeV. This latter pairing 
involves red (or color 3) light quarks left over above, and 
it further breaks the remaining color subgroup SU (2)c 
by selecting, say, scalar d^si, U3S2 pairs. What is finally 
left over in a "normal" state are S3 quarks: according 
to their perturbative pairing may still occur in a spin- 
1 channel pi^ . 

5. Summarizing this work, we first recall the lessons 
learned from Nc—2 QCD: here, diquarks are colorless 
baryons, and the scalar ones are as different from all other 
baryons as pions are different from other mesons. Chiral 
symmetry is restored at a very small chemical potential 

1 /2 

/i = 0{mq ), and the usual {qq) condensate is replaced 
by a diquark {qq) condensate. Our simple mean field ap- 
proach, lattice and instanton simulations are consistent 
with each other. 

In the real world with 7Vc=3, scalar ud-diquarks are 
strongly bound. If not for confinement, a Bose condensed 
diquark gas would be the ground state at low densities. 
The transition from nuclear matter to diquark-condensed 
matter starts at densities n;, ~ 0.7 — 1 fm~'^ , as a phase 
in which both (qq), (qq) are non-zero (it may, however, be 
unstable). An appreciable BCS-like gap A ~ 100 MeV 
builds up towards the chiral restoration point. Beyond, 
instantons become Debye-screened, so that a non-zero 
(qq) eventually only survives due to one-gluon exchange. 

An interesting way to understand these phenomena 
is suggested by comparing the features of the instanton 
ensemble at high T and high fi. In both cases quark 
propagation in time direction is favored over spacelike 
propagation (suppressed by exp(— TrTr) and exp(i/ir), re- 
spectively). As a result, the fermion determinant leads 
to strong correlations among instantons. Eventually, the 
random ensemble breaks into small clusters restoring chi- 
ral symmetry. These clusters are oriented in the time di- 
rection: "instanton-anti-instanton" pairs at high T and 
"polymers" at high density. 
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